Bosonic amplification of noise-induced suppression of phase diffusion 
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We study the effect of noise-induced dephasing on coUisional phase-diffusion in the two-site Bose- 
Hubbard model. Dephasing of the quasi-momentum modes may slow down phase-diffusion in the 
quantum Zeno limit. Remarkably, the degree of suppression is enhanced by a bosonic factor of order 
N/ log A'^ as the particle number A'^ increases. 
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The interplay between unitary evolution and decoher- 
ence has been a central issue of quantum mechanics. A 
universal formula has been put forward for the dynami- 
cal control of quantum systems weakly coupled to a bath 
[l[ . According to this formula, relaxation or decoherence 
can either be suppressed or enhanced by interventions 
whose rate is much higher than, or comparable to (re- 
spectively) the inverse non-Markovian memory time of 
the bath response. If these interventions are either pro- 
jective measurements or, equivalently, stochastic dephas- 
ing of the system's evolution , the resulting slow-down 
or speed-up of the relaxation/decoherence coincide with 
the quantum Zeno effect (QZE) [H 0] or the anti-Zeno 
effect (AZE) f?,!], respectively. 

This universal formula provides simple recipes for the 
dynamical projection of both single- and multi-partite 
quantum states, provided the bath spectral response is 
known. Yet an open question remains: can we similarly 
control/protect quantum states of large multi-partite sys- 
tems from the buildup of many body correlations among 
its N interacting particles ? The analysis of this scenario 
is tantalizing and nearly impossible for a multi-mode sys- 
tem with large N. However, useful insights can be gained 
by exploring few- mode models, for which full numerical 
solutions may be used to support analytic early-time ap- 
proximations. At such times, we may describe the slow- 
down or speedup of the system's many-body evolution 
by noisy perturbations, as QZE or AZE, respectively. 
Since phase-diffusion between atomic Bose-Einstein con- 
densates (BEG) 0, i, i, Q has a non-Markovian cor- 
relation time of ms, it is particularly amenable to the 
observation of the QZE and AZE 8]. Our main result 
here, is that for iV-boson condensates, the QZE is Bose 
amplified. 

Specifically, we consider the noise-induced suppression 
of phase-diffusion in the two-site Bose-Hubbard model, 
recently used to describe experiments of quantum inter- 
ference 0, Q and tunneling in an array of double wells 
[lo| . This model, under the tight-binding condition for 
3D wells, I ^ N\a\, where a is the scattering length, and 
I = ^Jh/imuQ) is the characteristic size of a trap with 
frequency wq, is accurately described by the quantized 




FIG. 1: (color online) Bose-enhanced QZE suppression of two- 
mode depletion: (a) Schematic diagram of the two-site Bose- 
Hubbard model. Depletion of the excited quasimomentum 
state is slowed down by a non-local noise source, (b) TV-partite 
Bloch sphere, with Lz-dependent depletion and L^-dependent 
noise, (c) Suppression of the depletion rate R as the overlap 
of Gdepiiijj) and Ft{uj) noise decreases with A'' (solid red to 
dashed red - Eqs. ([^ . (ITOll 1 . 



two- mode Josephson Hamiltonian here rewritten as 



H 



(1) 



where = {a[a2 + a\ai)/2, Ly — {a\a2 — alai)/{2i), 
and Lz — {fii — «.2)/2 generate the SU{2) Lie algebra. 
The mean-field values (Li) determine the reduced single- 
particle density matrix p^^^ — {a\aj)/N, with mode in- 
dices i,j = 1,2. The operators hi and a\ are bosonic 
annihilation and creation operators respectively, for par- 
ticles in mode i with corresponding particle number op- 
erators fii — a\ai. The bias potential is here set to 
zero, J is the intermode coupling, and U is the col- 
lisional interaction frequency. We have eliminated c- 
number terms, proportional to the conserved total parti- 
cle number N = fii -\- fi2. 
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The cohisional lI term in the Hamiltonian ([T]) leads to 
'phase-diffusion' @, 0] which degrades the reduced 

single-particle coherence. Its eigenstates, 
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|0), (2) 



constitute a preferred basis set resilient to this pro- 
cess [12]. On the other hand, the most sensitive states 
to phase-diffusion are the spin coherent states \6, (j)) = 
e~''^^^e~'^^!' |Z, — /), with = 7r/2, corresponding to equal 
populations of the two sites, and a well-defined relative 
phase (j). For fully separated modes (J = 0) the single 
particle coherence of these states is lost as exp [— (t/trf)^] 
with a characteristic decay time — {UVl)^^ and re- 
vives after U = tt/U 0, For finite J and U > 
(repulsive interaction) the fastest phase-diffusion occurs 
for the antisymmetric coherent state 
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\l,m) , (3) 



i.e. the state with all particles populating the excited, 
odd superposition of the modes. This experimentally 
realizable state 01 will be used as the initial condition 
throughout this work. Its evolution with [/ > is iden- 
tical to the evolution of the state |7r/2, 0) with [/ < 0, 
which for \UN\ > J drives the system towards a macro- 
scopic cat state Since [L ,H]= 0, we fix ^ = N/2. 

Control of the phase diffusion of the state Q may be 
attained by noise- induced dephasing of the odd- and even 
two- mode superpositions, as illustrated in Fig. [TJ This 
control can be introduced by any noise source that does 
not distinguish between the sites, e.g. an off-resonant 
incoherent light source focused at the barrier between 
them (Fig. [T^). Such noise, implemented perpendicu- 
lar to the depletion axis (Fig.[T}D), affects the stochas- 
tic (noisy) modulation of the splitting J, modifying the 
system Hamiltonian as 



Hs{t) = [J + hSs{t)]L, + ULl 



(4) 



The depletion-dephasing interplay is described by the 
general non-Markov (Zwanzig-Nakajima) second-order 
master equation (ME) for the reduced density matrix of 
the system. The resulting noise-controlled, normalized 
decoherence rate R conforms to the universal formula 
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expressing R{t) as the convolution of the finite time 
spectral intensity of the noise control (where et(w) = 



0.999 



0.998 
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FIG. 2: (color online) Single-particle coherence as a function 
of rescaled time, starting from the coherent state |7r/2,7r) 
with (a) K = 0.5, (b) it — 2. Bold lines in (a) correspond 
to iV = 100 (solid red), 150 (dashed blue), and 300 (dash- 
dotted green) particles. Bold lines in (b) are A'^ = 50 (dotted 
black), 100 (solid red), 200 (dashed blue), and 400 (dash- 
dotted green) particles. Gray lines correspond to the analytic 
forms of Eq. (0 for the same A^. 



dt'e'['^*'+-/'o''^=(*"' - )'**"l is the Fourier transform of 
the stochastic phase factor), and the Fourier trans- 
form of the depletion correlation function ^depi = 

(2/A^)2(e*'^^-(*-*')/'*L^(t)e-*"'^-(*~*')/'*i2(t'))- Equation 
([5]) provides a general recipe for controlling quantum de- 
pletion by noise. It yields the QZE limit of R{t) suppres- 
sion when Ft{Lu) is spectrally much broader than Gdepi 
(Fig. [TJ:), i.e. when the inverse width of these spec- 
tral functions (their memory times) satisfy {tc)noise 
{tc)depi- Conversely, it yields the AZE limit of R{t) 
enhancement when the two spectral widths or memory 
times are comparable and their spectral centers are mu- 
tually shifted. 

In what follows, we focus on the QZE limit, neglecting 
(tc)noise altogether, i.e. taking the broadband noise to 
be Markovian. In atomic BECs, this limit is obtained for 
l/(ic)noise > F^ > l/{tc)depi > kHz, being the rate 
of the Markovian dephasing. In this limit, we can use the 
Markovian quantum kinetic ME 
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Exact solutions of Eq. ^ can be found by expansion in 
the \l, m) basis set and numerical integration. In order to 
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analytically approximate the initial phase-diffusion, we 
truncate the hierarchy of dynamical equations for the 
Li operators, at second-order correlations to obtain the 
BogoHubov Backreaction (BBR) equations 0, for 
the mean- field single-particle Bloch vector s = 2{'L)/N 
and the correlation functions — 4:{{LiLj + LjLi) — 
2{Li) {Lj))/N'^ . Linearizing these equations around the 
|7r/2,7r) state, we obtain the initial dynamics of the nor- 

-1/2 



malized correlation function g[^l 

corresponding to the fringe visibility in interference ex- 
periments (21 • In the absence of noise [Fx = 0) we find, 

l-i^sin^(Ar) «<1 
[ai'Mr)) ={ , (7) 



Pl2 
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4cothf(2je) 2 
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sinh (At) k > 1 



where A — — k\ and tan 8 = \/l — k, k ~ UN/ J 
is the coupling parameter, and r = Jt is the rescaled 
time. Dynamical BEC depletion in the weak-interaction 
{k < 1) regime is thus bound and inversely proportional 
to the number of particles N. By contrast, for strong 
interactions (k > 1), the phase-diffusion rate is indepen- 
dent of the number of particles, but its onset time scales 
logarithmically with N . This strong interaction insta- 
bility may account for the rapid heating observed in the 
merging of two condensates with a tt relative-phase, on 
an atom chip [9]. 

Numerical results based on the Markovian ME confirm 
the short-time dynamics of Eq. ([7]). The weak-coupling 
behavior (Fig. [2^) exhibits the anticipated stable oscil- 
lations. The oscillation amplitude decreases with in- 
creasing N while keeping k fixed, so that 0{1/N) initial 
quantum fluctuations remain small compared to the 0(1) 
classical mean- field. By contrast, for strong interactions 
(Fig. ^p) , quantum fluctuations grow rapidly and single- 
particle coherence is lost. Equation ([7]) gives an accurate 
description of the depletion at short times and a good 
estimate for the phase-diffusion time. At longer times, 
depletion is significant and the linearized BBR equations 
are no longer adequate. The fully nonlinear BBR equa- 
tions, however, describe the loss of single-particle coher- 
ence with good accuracy. Partial and full revivals are 
observed due to the finite number {N + 1) of phase-space 
dimensions. 

We now proceed to explore the effect of noise on phase 
diffusion. It is evident from Eq. ^ that in the absence 




FIG. 3: (color online) Single-particle coherence as a func- 
tion of rescaled time, with k = 0.5 in (a),(b) and k — 2 
in (c),(d). The Hamiltonian dynamics with 7^ = (bold 
red lines) is compared to the evolution with non-vanishing 
7^ (normal blue lines). Solid, dashed, and dash-dotted blue 
lines in (a) and (c) correspond respectively to 7^, — 5, 10, 20 
and A'^ = 100. Normal blue lines in (b) and (d) corre- 
spond to 'yx = 10. Solid, dashed and dash-dotted lines 
in (b) correspond to = 100, 150, 300 particles, respec- 
tively, whereas solid, dashed and dash-dotted lines in (d) are 
= 100, 200, 400 particles, respectively. Gray symbols por- 
tray the QZE behavior of Eq. dS]). 



we can adiabatically eliminate Ayz in the linearized BBR 
equations and obtain the QZE behavior. 



QiMt) = exp ( -^t 



(9) 



where 'jx = ^x/J- The modified diffusion time in the 
presence of noise, = 2"fx/^^ = N^x tanh'^(2i8)/A^ ^ 
l/f2, should be compared, in the strong-interaction 
regime, with the noise-free, finite- J diffusion time, Td oc 
log [iVtanh^(2i8)] /A. Hence, the transition from the 
hyperbolic growth ([7]) which only depends on N through 
its onset time, to the QZE-suppressed depletion ([9]) at 
a rate linear in N, introduces a bosonic factor of order 
N/ logiV > 1 in the diffusion time ratio frf/r^. The QZE 
is thus strongly amplified with increasing particle num- 
ber. This constitutes the main result of this work. 

Similar scaling is obtained from the universal Eq. ([5]). 
In the QZE limit of non-Markovian depletion in the pres- 

(1)^ 



of noise, regardless of the interaction strength, g[^l at ence of Markovian noise, g\i{t) ^l-Jg R{t')dt' may be 
T < 1/A scales quadratically rather than linearly in time, inferred from it by substitution of the stochastic broad- 
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band Ft{uj) ~ (2/iV)2(t/27r)sinc^(cjt/2)(L2)^ go that 

2t f ijjt\ 

(8) R{t) - ^ / GdeA^ + J/h) sinc^ f - j {Ll)dw 



where = yj2/N\ coth(2i8)|A. When the frequent mea- 
surement dephasing condition 7^; 3> A > is satisfied, 



(10) 

This form characterizes frequent projective measure- 
ments H. The QZE is obtained when sinc^(wt/2)(L^) 
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FIG. 4: (color online) Absolute values of the elements of p 
in relative Fock basis \l,m), for = 30 and k = 2: (a) 
Initial coherent state, (b) cat state after noise-free evolution 
at T = 2.4, (c) incoherent mixture after evolution with local 
site noise with = 0.05J, at r = 2.3, (d) protected single- 
particle coherence after evolution with F^, = J. 

is much broader than Gdepi- As N is increased while 
keeping k fixed, we have 4:{L'^)/N'^ « ~ 1 and 
Gdepii^^ + J/K) o^l/N (Fig. [TJ;) so that R(t) cx 1/iV. 

The QZE suppressed phase diffusion is illustrated in 
the numerical results of Fig. [31 where we compare the ini- 
tial evolution oi g'^\ with and without noise, in the weak- 
and strong-interaction regimes. The weak-interaction os- 
cillations of Eq. d?!) are replaced, as "fx is increased, by 
the exponential decay of Eq. ([5]), at a rate proportional to 
l/{N"fx) (Figs. [3^, b). The strong-interaction dependence 
on the dephasing rate "fx (Fig. [3t) and its Bose-amphfied 
suppression (Fig. show a clear transition from log 
dependent diffusion-times followed by A^-independent de- 
pletion rate, to dependent depletion rates. These 
numerical results agree well with the appropriate closed 
form of Eq. ^ and Eq. ^. 

It is interesting to contrast this noise-induced QZE be- 
havior to the effects of local-site noise, which may be in- 
duced by collisions with thermal particles [l^. In Fig. [3] 
we plot the density matrix elements for the intial co- 
herent state (Fig. 2^) and for the macroscopic cat state 
generated after noise-free evolution in the strong inter- 
action regime (Fig. The dynamics leading to this 
cat state is extremely sensitive to local site noise [isj , 
which destroys the macroscopic coherence, resulting in a 
50-50 statistical mixture (Fig.^t). Phase-diffusion in this 
case, is enhanced to the extent that the strong-interaction 
diffusion-time is bound at large N By contrast, the 
site-indiscriminate (non-local) noise considered here pro- 
tects the single-particle coherence and slows down phase- 
diffusion (Fig. [Hi). Whereas any weak local noise will de- 
grade the intricate dynamics leading to a cat state, our 



non-local noise needs to be sufficiently strong to induce 
the QZE. 

To conclude, we have found novel collective features 
of the QZE, which do not appear in the noise-controlled 
decay of single particles. The results presented here are 
generic rather than specific to the two-site Bose-Hubbard 
model of atomic BECs. There is currently great inter- 
est in phase-diffusion experiments, enabling the measure- 
ment of single-particle coherence via the visibility of in- 
terference fringes 0, 0, 0|- Our predictions may thus 
be directly verified using current experimental apparata. 
Consequently, new avenues may be opened for noise- 
control of complex multipartite systems. 
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